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ABSTRACT 

We report the presence of two distinct red giant branches (RGBs) in the 
globular cluster NGC 288 from the narrow-band calcium and Stromgren h &c y 
photometry obtained at the CTIO 4m Blanco telescope. The RGB of NGC 288 
is clearly split into two in the hk [={Ca — h) — {h — y)] index, while the split is not 
shown in the h — y color. Unlike other globular clusters with multiple populations 
reported thus far, the horizontal branch of NGC 288 is only mildly extended. Our 
stellar population models show that this and the presence of two distinct RGBs 
in NGC 288 can be reproduced if slightly metal-rich {A[m/H] ^ 0.16) second 
generation stars are also enhanced in helium by small amount (Ay 0.03) and 
younger by ~ 1.5 Gyrs. The RGB split in hk index is most likely indicating that 
the second generation stars were affected by supernovae enrichment, together 
with the pollutions of lighter elements by intermediate-mass asymptotic giant 
branch stars or fast-rotating massive stars. In order to confirm this, however, 
spectroscopy of stars in the two distinct RGB groups is urgently required. 

Subject headings: globular clusters: individual (NGC 288) — stars: abundances 
— stars: evolution 



INTRODUCTION 



Unlike the conventional wisdom, observations made during the past decade have re- 



vealed that many globular clusters (GCs) are p ossessing more t 



tion. Some of these peculiar GCs, such as uj Gen (iLee et al. 



(iLavden fc Saraiedinilll997l : ISiegel et al.ll2007). M2 2 flDa Costa et al 



1999 



lan one stellar p opula- 
Bedin et al.ll2004[). M54 



2009 



Lee et al 



Marino et al.l l2009h . NGC 1851 faan et al.l l2009[ Ihee et al.l l2009b[ ICarretta et al. 



2009a 



2OIOI) 
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Terzan 5 flFerraro et alJl2009h . and NGC 2419 fICohen et al.ll20ld : lOi Criscienzo et al.1120111 ) 
show evidences of supernovae (SNe) enrichment, indicating that they are rehcs of more mas- 
sive primeval dwarf galaxies, rat her than beiiig nor mal GCs. For other GCs with mu ltiple 
pop ulations, such as NG C 2808 dPiotto et all boOTh . NGC 6388 dMoretti et allboogh . and 
M4 (IMarino et al.l 120081 ). however, the evidence for the discrete distribution of heavy ele- 
ments as observed in the RGB of u Cen is lacking, although spreads in some lighter elements 



(ICarretta et al. 



2005 



2009 



a nd references ther ein) and helium ( iD'Antona et al.l l2005l : iLee et al. 



Piotto et al.ll2007l : lYoon et al.ll2008l ) are reported. Therefore, the presence of chemical 



inhomogeneity and multiple populations in these GCs is largely considered due to the pollu- 
tion from the i ntermediate-mass asymptot i c giant branch (AGB) stars and (or) fast-rotating 
massive stars ( iVentura fc D'Antonall2008l : iDecressin et al.l 120071 ). which is expected even in 
normal GCs. 

Because of their important implications on the hierarchical merging paradigm of Galaxy 
formation, search for more GCs with dwarf galaxy origin (i.e., with evidence of SNe enrich- 
ment) would be extremely important. The purpose of this Letter is to report that NGC 288 
is also showing a clear split in the RGB from the narrow-band calcium photometry. This 
observation is compared with our stellar population models to argue that the two populatons 
are different in terms of overall metallicity, helium, and age by small amounts. 



2. OBSERVATIONS AND COLOR-MAGNITUDE DIAGRAMS 

Our observations in Ca, 6, and y passbands were performed using the CTIO 4m Blanco 
telescope on 2009 July 27. The telescope was equipped with eight 2048x4096 SITe CCDs, 
providing a plate scale of 0.27 arcsec pixel"^ and a field-of-view of 36x36 arcmin on the 
sky. However, since the 4x4 inch filters used in our photometry can not cover the entire 
field, only four CCD chips located in central region (i.e., chip 2, 3, 6, and 7) covering 
2048x2500 pixels per chip were used in the final photometry. The total exposure times for 
Ca, Stromgren 6, and y were 1650, 264, and 132 seconds, respectively, split into short and 
long exposures in each band. NGC 288 was placed on chip 6, approximately 3.0 arcmin 
South and 3.1 arcmin East from the CCD center. The IRAfQ MSCRED Package was used 
for preprocessing including bias correction and flat fielding. The brightnesses of objects in 
NGC 288 were measured with the point-spread function (PSF) fitting routine DAOPHOT 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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II and ALLFRA ME (IStetsonI Il987l . Il994l ). and aperture corrections were calculated using 
the DAOGROW Jstetsonlll99oR ~ Our photometry in Ca, b and y passbands were th e n use d 
to calculate the hk [= {Ca — b) — {b — y)] index defined by lAnthony-Twarog et al.l (Il99ll ). 
The Ca filter in the hk index is meant to measure essentially ionized calcium H and K 
lines, and the hk index is known to be about three times more sensiti ve to metallicity than 
^1 [= {v — b) — {b — y)] index is (ITwarog fc Anthony-Twarogl Il995l ). The same filter set 
empl oyed in this obs e rvation was extensivel y used by us in our previous investigations of 
GCs jRev et al.lboool . booi llee et aljl2009al lbll 

Figure 1 shows color-magnitude diagrams (CMDs) of NGC 288 in (6 — y, y) and {hk, 
y) planes. To examii ie the CMP featur es more clearly, magnitude error, chi, sharpness 
and separation index (IStetson et al.ll2003l ) were used to reject stars with large photometric 
uncertainty and those affected by blending and adjacent starlight contamination. All stars in 
the Figure 1 lie within the chip 6, and therefore our CMDs are not subject to any uncertainty 
stemming from the possible chip to chip variations of the mosaic CCDs. The most remarkable 
feature of Figure 1 is the presence of two distinct RGBs in the hk vs. y CMD. When measured 
aX y = 16.5 mag, the mean separation between the two RGBs is about 0.10 mag in hk index. 
The discrete distribution shown in RGB, however, is not apparent in the subgiant branch 
(SGB). Note also that the RGB split is not shown in {b — y, y) CMD. This is most likely 
because the Ca filter in the hk index is much more sensitive to changes in Ca abundance 
than other color indices like b — y. 

Given the small foreground reddening value of E{B — V) = 0.03 ( HarrisI 1996 ) to- 
ward NGC 288, it is very unlikely that the differential reddening has caused the double 
RGBs. Furthermore, in contrast to other color indices, the hk index is known to be in- 
sensitive to interstellar reddei ng, E{hk)/E{B — V) = —0.12 and E{hk)/E{b — y) = —0.16 
(lAnthony-Twarog et al.lll99ll ). Therefore, if we adopt the reddening of E{B — V)= 0.03, 
E{hk) = —0.0036 is obtained for NGC 288, which is negligible compared to the separation 
in hk index (~ 0.10 mag) between the two RGBs. Star counts of two subpopulations indi- 
cate that the bluer RGB population ("Pop-1") takes about 60% of the whole population, 
while the redder RGB population ("Pop-2") comprises about 40% of total population, in the 
magnitude interval y = 16.5 ±1.5 mag. This ratio is not sensitive to the adopted separation 
index in our photometry. 



3. COMPARISON WITH STELLAR POPULATION MODELS 



In order to better understand the origin of the RGB split in hk index, and to place 
constraints on the chemical combinations of two subpopulations, we have constructed stellar 
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popu lation models based on the latest version of the Yonsei-Yale (Y^) isochrones (lYi et al. 
20081) and HB evolutionar y tra cks (S.-I. Ha n et al. 2011, in preparation). Readers are referred 



to iLee et al.l ( Il990l . Il994j ) and lYoon et al.l ( 120081 ) for the details of our model construction. 
Figures 2 and 3 present our synthetic CMDs for NGC 288 in {hk, y) and {b — y, y) planes, 
respectively. Our models were constructed under three different assumptions regarding the 
chemical enrichment and age spread in NGC 288. First, we assumed that the second genera- 
tion population (Pop-2; redder RGB) is more enhanced in metallicity and younger (A[m/if] 
^ 0.16 dex. At ^ 1.^ Gyr), but not enhanced in helium abundance [hereafter AZ+AAge 
model; panel (b) in Figures 2 and 3]. These models match well with the observed CMDs from 
the MS through the RGB in {hk, y) and (6 — y, y) planes. Yet, the models fail to reproduce 
the HB, as the synthetic HBs are too extended in color including significant numbers of RR 
Lyraes and red HB stars (HB typ^ = 0.61), while the observed HB is only mildly extended 
with mostly blue HB stars (HB type = 0.91). This is because both metal enhancement and 



younger age in second population move the HB to red in CMD (see lLee et al.lll994j ). Second, 
we then assumed that both metal and helium abundances are enhanced (A[m/if] ~ 0.16 
dex. Ay ^ 0.03 dex), but age is constant [hereafter AZ+AY model; panel (c) in Figures 2 
and 3]. These models match well with the observed CMDs from the RGB through the HB 
in (/ifc, y) and (6 — y, y) planes. However, they can not reproduce the narrow SGB in {hk, 
y) CMD. Therefore, both AZ+AAge and AZ+AY models are in conflict with the observed 
CMDs of NGC 288. 

Finally, we assumed that not only metal and helium abundances are enhanced, but also 
age is younger in Pop-2 [hereafter AZ+AY+AAge model; panel (d) in Figures 2 and 3]. 
These models are in good agreements with the observations from the MS to the HB. The 
enhanced metal abundance in Pop-2 makes the RGB split in hk index as observed, while the 
younger age can explain the narrow and apparentl y single SGB in (hk, y) CMD. The increase 
in helium abundance in Pop-2 moves HB bluer (see lLee et al.ll2005l ). almost cancelling out the 
effects by enhanced metallicity and younger age, making the blue HB only mildly extended 
with the HB type similar to the obs erved value (HB type = 0.90). Note that, in our HB 
simulations, we employ the standard iReimersI ( 119771 ) mass-loss law and the same mass-loss 
parameter t] for the two subpopulations. Input parameters adopted in our best models (i.e., 
AZ+AY+AAge model) are listed in Table 1. Our models are computed with the same 
abundance of [CNONa/Fe] for both Pop-1 and Pop-2. More detailed models including the 
possible difference i n [CNONa/Fe] bet ween the two subpopulations would change the age 



estimates (see, e.g.. ICassisi et al.ll2008l ). 



^The HB type is the quantity, (B-R)/(B+ y+R). where B V, and R are the numbers of bhie HB, RR 
Lyrae variable, and red HB stars, respectively ([Lee et al.lll994l ). 
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4. DISCUSSION 

We have shown that the RGB of NGC 288 is spht into two distinct sequences. While 
this is most hkely the effect of Ca II H & K lines, it is important to check whether the CN 
band at A ~ 3870 A could affect the hk index due to the proximity o f the CN band to th e 



blue tail of Ca filter transmission curve. In particular, according to iKayser et al.l ( 120081 ). 
lower RGB stars in NGC 288 exhibit CN bimordality that spans about 0.6 dex. In Figure 
4a, we have matched their spectroscopic data with our photometry, where we can see that 
"CN-stong" stars lie well on the redder RGB sequence, whereas "CN-weak" stars are on 
the bluer RGB. This further suggests that possible contamination of CN band to hk index 
should be investigated in more detail. It is known that the passband of the narrow band 



Clarke et al. 


1975; 


Lee et al. 


2009a) 



of incidence beam (see 
relevant when the Ca filter is used with a relatively fast telescope like the prime focus of the 
CTIO 4m telescope, where the central wavelength drift of the Ca filter is estimated to be 
about 15 A to th; shorter wavelengtlfi (see Figure 4b). 

In order to see the influence of CN band on the hk index in our photometry, we have 
calculated synthetic spectra using the ATLAS 9 model atmosphere (ICastelli fc Kuruczll2003l ) 
for the RGB star at the magnitude level of HB with Teg = 4750 K, logg = 2.0 (in cgs unit), 
■^turb = 2.0 km/s, and [Fe/H] = —1.6 (see Lee et al. 2009a Supplementary Information 
for detail). The typical RGB stars in GCs show an anticorrelation between CN band and 
CH band strengths and a correlation between C N band and NH Ijand strengths, indicating 
that the nitrogen controls the CN band strength (IBriley &: Smithlll993l ). Therefore, we have 
compared the synthetic spectra between the stars with normal and enhanced nitrogen abun- 
dances (see Figure 4c). We obtain Ahk ~ +0.006 for A [N/Fe] = +1.0 dex, which suggests 
that the influence of CN band on the hk index would be negligible in our photometrj]^. 
Recently, ISbordone et al.l (120111 ) also reached at a similar conclusion, finding that CNONa 
anticorrelations have an effect of at most ~ 0.04 mag in the hk color at fixed M^. Their 
calculations include not only the direct effects of CN and other absorption bands (NH & 
CH), but also the effects on the continuum levels in the Ca, b, and y filters. 

If the CNONa abundances have only little effects on the hk index as discussed above, 
which should be confirmed in the forthcoming works, the RGB split discovered in our hk vs. 
y CMD would indicate a small difference in Ca abundance between the two subpopulations. 



http://www.ctio. noao.edu/instruments/filtcrs/filtcrs_66. html 

^If we have assumed the wavelength drift of 20 A, a similarly small value of Ahk (^ +0.009) is obtained, 
confirming that this result is not very sensitive to the uncertainty in the value of the wavelength drift. 
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Since calcium and other heavy elements can only be supplied through SNe explosions, this 
in turn would suggest that the second generation stars were affected by SNe enrichment, 
together with the pollutions of lighter elements (such as the enhanccnmcnt of N and the 
depletion of O) by intermediate-mass asymptotic giant branch stars or fast-rotating massive 
stars. Spectroscopy of stars in the two distinct RGB sequences is crucial to confirm the small 
difference in the abundance of heavier elements suggested in our photometry. 

Wc thank the referee for a number of helpful suggestions. Support for this work was 
provided by the National Research Foundation of Korea to the Center for Galaxy Evolution 
Research. This material is based upon work supported by AURA through the NSF under 
AURA Cooperative Agreement AST 0132798, as amended. 
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Fig. 1. — Color-magnitude diagrams for NGC 288 from CTIO 4m Blanco Telescope. Note 
the discrete double RGBs in the hk vs. y CMD. Two RGBs are separated by ~ 0.10 mag in 
hk index. The reddening vector values are too small [E[hk) — —0.0036, E(b — y) — 0.0225, 
and Ay — 0.093] to present in this figure. 
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Fig. 2. — Our population models for NGC 288 compared with the observed CMD in hk vs. 
y plane. In the models, the red lines and symbols are for the metal-rich second generation 
population, while the blue lines and symbols are for the metal-poor poupulation. The crosses 
in the models denote RR Lyrae variables. The best match is obtained when metal-rich 
{A[m/H] f« 0.16) second generation stars are also enhanced in helium by small amount (Ay 
~ 0.03) and younger by ~ 1.5 Gyrs (see text). 
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Fig. 4. — (a) Same as Fig. lb, but zoomed around the SGB and RGB regions in hk vs. y 
CM P. Denoted by blue a nd red dots are respectively "CN normal" and "CN strong" stars 
from iKayser et al.l (120081 ). The CN strength seems to trace the calcium abundance, (b) A 
comparison of Ca filter transmission functions. The Ca filter at the prime focus of the CTIO 
4m telescope is shifted approximately 15 A to the shorter wavelength compared to the case 
of perfectly collimated beam, (c) Synthetic spectra for the CN normal (the blue line) and th e 
CN strong (the red line) RGB stars (see also Supplementary Fig. lb of iLee et al.ll2009al ). 
Although the CN band starting at A ~ 3885 A hes on the blue tail of Ca filter transmission 
function, contamination from the CN band appears to be negligible (see text). 
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Table 1. Input Parameters Adopted in Our Best Simulation of NGC 288 



Parameter 


Population 1 


Population 2 


Z 


0.00083 


0.00121 


Y 


0.231 


0.258 


[a/Fe] 


0.3 


0.3 


Age 


13.7±0.3 Gyr 


12.2±0.3 Gyr 




0.53 


0.53 




0.2179 


0.2094 




0.020 


0.020 


Population Ratio 


0.6 


0.4 



Reimera ( 119771 ) mass-loss parameter. 
^Mean mass-loss on the RGB (Mq). 
'^Mass dispersion on the HB (Mq). 



